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We present a theoretical study of the resonant interaction between dynamical localized states 
{discrete breathers) and linear electromagnetic excitations (EEs) in Josephson junction ladders. By 
making use of direct numerical simulations we find that such an interaction manifests itself by 
resonant steps and various sharp switchings (voltage jumps) in the current-voltage characteristics. 
Moreover, the power of ac oscillations away from the breather center (the breather tail) displays 
singularities as the externally applied dc bias decreases. All these features can be mapped to the 
spectrum of EEs that has been derived analytically and numerically. Using an improved analysis of 
the breather tail, a spectroscopy of the EEs is developed. The nature of breather instability driven 
by localized EEs is established. 
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I. INTRODUCTION 

Various nonlinear and discrete systems (nonlinear lat- 
tices) have attracted a lot of interest as they display di- 
verse fascinating phenomena.EJ Well known examples of 
such phenomena are solitary excitations, propagation of 
(non)linear waves and appearance of various inhomoge- 
neous structures. 

Moreover, the interest in thisparea was boosted by the 
prediction, jlheoretical analysi^Q and the subsequent 
observationErtj of intrinsic dynamic localized excitations 
(discrete breathers) that are periodic in time and local- 
ized in space. Note here, that the origin of such dynam- 
ical localization is not the presence of disorder but the 
interplay between the nonlinearity and discreteness. 

These peculiar states have been experimentally werified 
as vibrational modes in low-diniepsional crystals,El local- 
ized excitations in spin latticesjj aH|d-4ocalized resistive 
states in Josephson junction arraysLnld. The latter sys- 
tems are of special interest because they have served for 
many years as well-controlled l|abDratory objects to study 
various nonlinear phenomena.EHlj Moreover, at variance 
with other systems intrinsic localized modes found in 
Josephson coupled systems can be excited in the pres- 
ence of time independent external driving forces. 

A well known structure where dynamical localized 
states appeasj, is the anisotropic Josephson junction lad- 
der (JJL)Lnl2l. A schematic view of such a ladder is given 
in Fig.|l]. The ladder contains small Josephson junctions 
indicated by crosses in Fig.^ in both longitudinal [ver- 
tical junctions) and transverse (horizontal junctions) di- 
rections to the dc bias current 7. The anisotropy of the 



ladder is due to the different sizes of vertical and hori- 
zontal junctions and is characterized by the anisotropy 
parameter ry = j^, where Icv and Ich are respectively 
the critical currents of the vertical and horizontal junc- 
tions. 
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FIG. 1. Josephson junction ladder. Crosses mark the 
individual junctions. Arrows indicate the direction of external 
current flow (dc bias 7). 

Dynamical localized excitations persist in a JJL due to 
the intrinsic bistability property of a single small under- 
damped Josephson junction. One of these stable states is 
a superconducting one with zero voltage drop across the 
junction. The other state is a resistive one with a nonzero 
voltage drop, also called whirling state. A breather state 
in a JJL is characterized by a few junctions being in the 
resistive state, while the rest of all junctions are in the su- 
perconducting state. The presence of breather states can 
be verified by measurements of a total dc voltage drop 
across the ladder, which is used to plot current-voltage 
(/ — V) characteristics. This method does not provide 
spatially resolved information. It has been successfully 
combined with snapshots made using low temperature 
laser microscopy techniquesj3 which allow for a spatial 
resolution of the dc voltage drops. Note that both mcth- 
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ods provide only with time-averaged voltage drop data, 
so the internal dynamics are so far not accessible in ex- 
periments. 

However, the full dynamical picture is much more sub- 
tle than the time- averaged picture might suggest. In par- 
ticular the Josephson junctions in the superconducting 
state exhibit small librations of the Josephson phase and 
correspondingly, nonzero ac voltage drops. The ampli- 
tude of these librations should decay to zero with increas- 
ing distance from the resistive junctions. It is also well 
known that various systems of coupled Josephson junc- 
tions support a delocalized class of excitations, panicly, 
small amplitude electromagnetic waves (EWs) .E2JI13 In 
J JLs the spectrum of EWs consists of three branches and 
depends in a complex manner on the anisotropy rj and the 
dc bias 7. The resonant interaction of these EWs with 
the homogeneous whirling state (HWS) in the presence 
of an externally applied magnetic field has been studied 
in [Ref. |l^. It was shown that such an interaction leads 
to resonant steps in I ~ V curves and the voltage posi- 
tions of resonant steps can be mapped onto the spectrum 
of EWs. |-| 

Early theoretical studiesll^ dealt with the possibility 
of resonant interaction of breather states with EWs of 
the ladder. Due to the intrinsic spatial inhomogeneity 
of a breather state, EWs can be excited even in the ab- 
sence of an externally applied magnetic field. The res- 
onant interaction of EWs with the breather state mafe 
ifests itself through the appearance of resonant stepst2l 
and various switchings between different breather states 
{"voltage jumps") in I — V curves. Moreover, the ampli- 
tude of the Josephson phase librations at some distance 
from the breather center increases drastically due to this 
interaction.E^I 

The spatial inhomogeneity of a breather state allows 
also for the appearance of localized small amplitude elec- 
tromagnetic excitations (EEs). We will show that the 
presence of localized EEs and their resonant interaction 
with the breather is of crucial importance. Our study 
resolves a long standing puzzle of the nature of breather 
instability. We show that most of these instabilities are 
driven by localized EEs. Especially the so-called retrap- 
ping (i.e. the switching from a breather to the supercon- 
ducting state) is due to these instabilities and can not be 
explained using standard retrapping arguments. 

In this paper we present a consistent theoretical study 
of resonant interactions between the breather states and 
EEs. We will derive the spectrum of EWs and calculate 
the ac power of oscillations at some distance from the 
breather center ("breather tail"). By making use of di- 
rect numerical simulations of the dynamics of JJLs in dif- 
ferent parameters ranges, we demonstrate that the I —V 
curves display a variety of different resonance steps and 
voltage jumps as the external dc bias, and correspond- 
ingly the breather frequency, decrease. All these features 
are mapped onto various resonances of the breather fre- 
quency (or even its second and third harmonics) with EE 
frequencies, as well as with combinations of frequencies 



of different EEs. Moreover, by monitoring the power of 
ac oscillations as a function of the dc bias 7, we develop 
a spectroscopy of the EWs in the JJL. 

The paper is organized as follows: In section II we 
derive the equations of motion within the frapiework of 
the Rcsistively Shunted Junction (RSJ) modelEiJ and ob- 
tain the spectrum of linear EWs of a JJL. A symmetry 
classification of different types of breathers is presented 
in section III. An improved analysis of the breather tail 
and the corresponding dependence of the power of ac os- 
cillations on the breather frequency is given in section 
IV. In section V we will classify different resonances of 
breathers with EEs. Direct numerical simulations of the 
dynamics of breather states are presented and discussed 
in section VI. Section VII is devoted to the interpretation 
of the obtained resonances and switchings. 

II. DYNAMICS OF A JJL AND THE SPECTRUM 
OF LINEAR ELECTROMAGNETIC WAVES 

The complete dynamics of a JJL is determined by the 
time dependent Josephson phases of vertical c/)^, upper 
horizontal (j)'^ and lower horizontal ij)!^ junctions. The 
subscript n labels the cell nurpber. By making use of the 
RSJ model for each junctionllil we obtain the following 
set of equations: 

mK) = i: 

m<Pn) = \ln (1) 

where the nonlinear operator Af is defined as 

Miy)=^y + ay + siny . (2) 

Here, the unit of time is the inverse plasma frequency 
Wp"^, and the currents I'l and are measured in units 
of the critical current of vertical junctions. The dimen- 
sionless parameter a determines the damping strength in 
each junction. Note that the positive current direction is 
chosen to be directed from bottom to top and from left to 
right. The currents flowing via the Josephson junctions 
and the Josephson phases are governed by the Kirchhoff 
laws 

TV , jh jh 

! 71 ' 71 71 — 1 

^ = r::-iil + iiU, (3) 

and the flux quantization law in each cell: 

- = € + C+1 - -^^ - C • (4) 

Here, we introduced the mesh currents /™ and the nor- 
malized inductance of the cell, (3l ■ 

The Kirchhoff equations can be subtracted from each 
other yielding 
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c 



(5) 



where C is a constant for the whole ladder. This constant 
corresponds to the net difference between the currents 
flowing through the upper and lower horizontal junctions. 
For an open ladder of finite size, C is zero. For a ladder 
of annular geometry with periodic boundary conditions 
C may be nonzero and corresponds to the flux " trapped" 
by the ladder ring. In the following we will consider the 
case of a finite open ladder with C = 0. Then we may 
always eliminate the currents from the set of equations 
(|l|) as — —In- Since the jimction width is larger tha-Q 
the London penetration depth, the mesh currents arellil 

C = In , (6) 

and the currents flowing through the vertical junctions 
are expressed as 



JV 



7- 



(7) 



Inserting the relations (Q), (^, and (0) into we 
finally obtain— jthe following set of coupled differential 
equationaa'tHlEjIlj : 

C + a^n + sin^r. = 7 + -^(AC + V0^i - V4>n-i) 

PL 



VPL 

C + + sin = — (V(/);; + cj)^-(j)^), 



I 



(8) 



where we use the notations A/„ = /„-i — 2/„ + /„+i and 

V/n = fn+l — fn- 

Next, we carry out the analysis of the delocalized class 
of excitations, namely of small amplitude electromagnetic 
waves (EWs). Note here that the spectrum of EWs oj{q) 
depends crucially on the state of the system. In the fol- 
lowing we will consider a static (superconducting) state, 
0*" — arcsin7 and 0*'* — 0*'' — 0. We decompose the 
Josephson phases into the particular form: 



(9) 



where ip"^, tf^ and (^J^ describe the small amplitude EWs. 
Substituting these expressions into system (^) and using 
the smallness of the amplitude of EWs we obtain 

+ a^l + s/l^^l = -^(A< + V^ti - V<^^i) 

PL 

(pi + oc^t + = -^{'^vl + - ^l) (10) 

WL 

In the weakly damped case as the parameter a 1 we 
can derive the spectrum uj{q) of EWs neglecting effects 



of dissipation. By taking the Josephson phases (p^, (/sjj 
and (p^ in the form 





(11) 



we find that the spectrum consists of three branches. The 
first is given by 



= 1, A,„ = 0, A,, = A,, 



(12) 



This branch is dispersionless and EWs corresponding to 
this branch are characterized by nonactive vertical junc- 
tions and in phase (symmetric) librations of the Joseph- 
son phases of upper and lower horizontal junctions. 

The two other solutions are generalizations of those 
discussed in [Ref. 13|, namely 



^F± y/F^ - G 



F 



1 
2 



+ — + ov/l-7' + -^(l-cosg) , (13) 

. PL 



G ^ (1 + — )v/l - 7' + -^(1 - cosg) . 



PlV 



Pl 



Both branches have a nonzero dispersion. 

The branch u;+ is characterized by A/j = — A/j for all 
wave vectors g, i.e. the upper and lower horizontal phases 
are antisymmetric. The frequency range of the branch is 
above the degenerate branch wq, i-e. w+(g) > ujq and it 
depends strongly on As the parameter Pl increases, 
the width of w+ (g) decreases and the branch approaches 
the dispersionless one, wq. In the opposite case of small 
(3l, the frequencies ijJ+(g) increase as l/V^. For zero 
wave number g = 0, the amplitudes of EWs in this branch 
are characterized by A.;, = and A^ — — A^, which 
means that only horizontal junctions are excited. 

The branch uj- becomes dispersionless for the particu- 
lar case of 7 = and it corresponds to the dispersionless 
band obtained in [Ref. |3). The frequency range of this 
branch is located below wo i-c ^-{q) < ^oil)- For zero 
wave number q — the horizontal junctions are not active 
{Ah — Ah — 0) and only vertical junctions are excited. 

For a finite size ladder with open boundary conditions 
and N cells, i.e. N + 1 vertical junctions, the spectrum of 
linear waves is discrete and characterized by the following 
choice of allowed wave number values: 



It: 



qi 



N + 1 



Z = 0,1, 2, ...,7V. 



(14) 



These EWs are cavity modes of the JJL. Odd values of I 
correspond to antisymmetric eigenvectors (with respect 
to reflections at the center of the ladder), whereas even 
values correspond to symmetric ones. 

Note that the above spectrum of EWs ( |l3|) is in gen- 
eral quite different frorp_.the EW spectrum of the homo- 
geneous whirling statell3 The latter can be obtained by 
choosing 7 = 1 in (Oh. The main difference is the ap- 
pearance of a gapless(acoustic) lower branch. 
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III. SYMMETRIES AND DC BIAS DEPENDENT 
FREQUENCIES OF BREATHERS 

In this Section we turn to the analysis of dynamic lo- 
calized excitations (breathers) in JJLs. As mentioned 
in the introduction, breather states correspond to a few 
junctions being in the resistive state with aU other junc- 
tions being in the superconducting state. The Josephson 
phases of resistive junctions are unbounded in time and 
the Josephson phases of superconducting junctions dis- 
play small amplitude librations with a frequency f2. This 
frequency is ca|llad,the breather frequency. 

ExperimentsQ ciJ have revealed many different breather 
structures. All of them can be classified into three sym- 
metry types using the reflection symmetries of the JJL. 
Some possible realizations are presented schematically in 
Fig. 2. 



a) b) c) d) 

FIG. 2. Examples of different types of breathers: a) 
up-down symmetry, b) left-right symmetry, c) inversion sym- 
metry, d) no symmetry. Black spots indicate the positions of 
whirling junctions. 

Breathers from the first group reveal an 'up-dovifn' 
reflection symmetry Sud (see, Fig.||a), i.e. they are 
invariant under exchange of upper and lower horizontal 
junctions. The second group consists of breathers invari- 
ant under a 'left-right' reflection symmetry Sir (see, 
Fig.^), i.e. they are invariant under a reflection at a 
vertical line cutting the ladder (this line is located either 
in the middle between two vertical junctions, or passes 
directly through one vertical junction). A third distinct 
^roup of breathers possesses an 'inversion' symmetry 
Sin (Figj2[;), i.e. these breathers are invariant under a 
reflection at a point which is either located on a vertical 
junction or in the center of a plaquette. A fourth group 
of breathers has no symmetries at all and does not 
belong to any of the three listed symmetry types. A par- 
ticular example of a breather without symmetry is shown 
in Fig||d. All of these types of breather excitationa have 
been observed experimentally and numerically.Lrli3 Each 
group of breathers can also have a different number of 
vertical junctions in the resistive state. Note that the 
particular example in Fig.||a possesses not only Sud sym- 
metry, but also Sir and Sin symmetries. However it is 
also possible to construct more complex breather states 
which display Sud symmetry only. 

Next we derive the average voltage drop across a re- 
sistive vertical junction V = (</>") for different breather 
types. For the particular case of a breather with up-down 
symmetry Sud, there are k rotating vertical junctions in 
cells (i -I- l)..(j -f- k) and two rotating horizontal junctions 



in the (i)th and (i -I- fc)th cells respectively. By making 
use of the flux quantization law (|4|) we obtain identical 
voltage drops across the vertical junctions in the resistive 
state. For the same reasons the voltage drops across the 
horizontal junctions are two times smaller. Neglecting 
nonlinear contributions from the time average of sin cf) on 
resistive junctions we obtain: 



aV = J + {III,) - {ID 
aV = J + {III,) - {I 



(15) 



i+fe-l/ 



Thus, the voltage drop across a resistive vertical junction 
that corresponds to the experimentally measured voltage 
drop across the ladder, is given by 



V = 



k"/ 



a{k + ri) 



(16) 



Similarly we analyze a breather with left-right symme- 
try Sir and with inversion symmetry Sin- Taking into 
account that in these cases the voltage drops across re- 
sistive horizontal and vertical junctions are identical, we 
find 



V 



kj 



i{k + 2ri) 



(17) 



In a similar manner the result for a breather which has 
no symmetry (cf. Fig. 2d) reads 



V = 



kj 



l(fc+§77) 



(18) 



The above results for the dependence of the average volt- 
age drop on the dc bias may be combined in a single 
expression 



V = 



k"f 



a[k + {3 - ^S)r)] 



(19) 



where k is the number of vertical rotating junctions and 
6 denotes the number of resistive horizontal junctions. 
Note that (5 = 4 for breathers with up-down symmetry, 
6 — 2 for left-right or inversion symmetry, and 6 = 3 for 
no symmetry. 

In order to analyze the interaction between the 
breather state and the linear EEs we need to know the 
frequency SI of a breather solution. Noting that the 
breather frequency is given by the lowest realized voltage 
drop across a resistive junction, we find 

1) up-down symmetry 



kj 



2a{k -I- ri) 



(20) 



2) left-right symmetry and inversion symmetry 



4 



n = 



a{k + 2ri) 



3) no symmetry 

n 



a{2k + 377) 



(21) 



(22) 



IV. SPATIAL TAILS OF BREATHERS 

We consider the spatial dependence of Josephson 
phases in the presence of a breather state. At some 
distance from the breather center ("breather tail") the 
Josephson phases librate with small amplitudes. In or- 
der to analyze the breather dynamics in the tail, we use 
the linearized system of equations (^0|) . Keeping in mind 
the time periodicity of the breather solution, the librat- 
ing Josephson phases take the form (for n < 0, and with 
the breather center as the origin) 




Xn-\-int 




(23) 



Introducing new sum and difference variables of the am- 
plitudes of Josephson phases of horizontal junctions 



with 



A = ln(z + - 1) 



4 - PlrjAB 
4 - 2f3LVB 



(28) 



(29) 



Note here that this expression can be obtained directly 
from ( p^ ) by assuming that q = iX and substituting Q'^ — 
iail instead of lu^. . 

The real 5R(A) and imaginary 5(A) parts of A determine 
respectively the spatial decay and spatial period of oscil- 
lations of Josephson phases in the breather tail. More- 
over, 3?(A) and 3(A) strongly depend on the breather 
frequency that in turn, can be changed by varying the 
external dc bias 7. In Fig||(a,b) we plot the real and 
imaginary parts of A for three breather types (cf. Fig. |^) 
versus 7. The minima of the real part of A correspond to 
resonances with linear EWs uj+ (q) and co- (q). As one can 
see, the current positions of these minima shift for differ- 
ent types of breathers. Note that the 7 range in the plots 
extends from zero to one. In fact the breather will ex- 
ist only in a narrower current region. This is due to-the 
presence of both finite nonzero retrapping currentli2l'E£l 
(switching to the superconducting state) and a partic- 
ular current 7 < 1, where the breather switches to the 
HWS. 



A^ = ^(A. 
it follows that 



A/.) 



(24) 



2 2 

(A - — coshA)A„ - -5-(l - e" 

PL PL 
PLV 



)A,; = 



SA,; = 



(25) 



{-n^ + ian + i)A+ = , 

where the frequency dependent parameters A and B are 

2 



Pl 



2 



From 

Ar = 



Sp we immediately find that A^ = 
\.h- It follows that ip'^ = —(pn, 



(26) 

and hence 
.e. breather 



tails appear with perfect up-down symmetry. This is at 
variance with the complex symmetry properties of the re- 
sistive breather center (see the Section III) . A nontrivial 
solution to the first two equations in (25) exists if 



AB 



Pl 



-B cosh A 



cosh A = 



(27) 



The dependence of the complex parameter A on the 
breather frequency Q is given by 
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FIG. 3. Dependence of (a) real part 5R(A), (b) imaginary 
part 9(A), and (c) Pac on 7 for different types of breathers: 
solid line - up-down symmetry, dashed line - left-right symme- 
try, dotted line - no symmetry. The parameters are a = 0.1, 
fiL = 3.0, 77 = 0.35 and fc = 1. 
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Since the EW frequencies uj+ {q) decrease with increas- 
ing inductance of the cell (|l^) , the position of the global 
minimum also depends strongly on 0r. (see FigHfa)). 
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FIG. 4. (a) K(A) and (b) P^c versus 7 for different values 
of the inductance of the cell: I3l = 3, solid line; /3l ~ 1, 
dashed line; /3_l = 0.2, dotted line. All of these curves are for 
a left-right symmetry breather with k — 1. 

Although the dynamics of the Josephson phases in the 
breather tail are completely determined by the depen- 
dence of the parameter A on the breather frequency, it 
may be more convenient to measure just the time-average 
power Pac of the libration of a junction at the edge of a 
JJL. We show below that the monitoring of Pac upon de- 
creasing the dc bias allows one to develop a spectroscopy 
of EWs. The value of Pac is determined by the average 
kinetic energy of the edge vertical junction 



(30) 



and is derived in Appendix A. The typical dependences 
of Pac on the dc bias and the inductance of the cell for 
different types of breathers are presented in Figs. ||c;, 



V. BREATHER INSTABILITIES AND 
CLASSIFICATION OF RESONANCES 

In early experiments on breathers in JJLsl'i it was ob- 
served that breather states may switch to other breather 
states or HWSs upon variation of the dc bias 7. These 
switchings can be either due to the disappearance of the 
breather state as a solution of the underlying dynamical 
equations or due to the novel effect of dynamical instabil- 
ity of the breather state. In the latter case the solution 
continues to exist, but turns unstable, forcing the sys- 
tem to search for another stable attractor state. The 



cause of such instabilities is the resonant interaction of 
the breather state with small amplitude (localized or de- 
localized) EEs. 

In order to analyze the breather instabilities in more 
precise terms, we have to linearize the phase space flow 
around a given breather solution and study the corre- 
sponding Floquet eigenvalue problem (see Appendix B). 
The outcome of this analysis is a spectrum of Floquet 
multipliers (i.e. eigenvalues) fi. In most cases the multi- 
pliers reside on a circle of radius 



(31) 



in the complex plane, which is less than 1. Thus, the 
multipliers can be expressed in the general form 



fi — R(a) exp ztiuj 



2n 



(32) 



where n is the breather frequency and uj is some char- 
acteristic EE frequency. The corresponding eigenvectors 
may be divided into two classes, namely those which are 
localized on the breather and those which are delocalized. 
While it is more involved to make analytical predictions 
for the class of localized Floquet eigenstates, we may im- 
mediately proceed with the characterization of the class 
of delocalized eigenstates. Although the concrete form of 
the delocalized eigenvectors has to be obtained numer- 
ically, the fact that they are delocalized allows one to 
determine their frequency to. Indeed, since the breather 
is itself a localized state, delocalized Floquet eigenstates 
simply correspond to the above discussed linear EWs. 
Their frequencies have been derived in Section II, and 
thus we can reconstruct the delocalized part of the Flo- 
quet spectrum using them. 

Stable breathers are characterized by all Floquet mul- 
tipliers being located inside or on the unit circle. Insta- 
bilities occur after collisions of multipliers on the inner 
circle with radius ( |3l] ) and a subsequent detaching from 
this circle towards larger absolute values. Although in 
the dissipative case an additional change of the control 
parameter (dc bias 7) is necessary in order for the cor- 
responding Floquet eigenvalue (s) to cross and escape the 
unit circle, we may still look for collision conditions and 
classify possible instability scenarios. 

Note that for Hamiltonian systems any collision and 
detaching from the unit circle leads to an instability. Fi- 
nite (even though possibly weak) dissipation may dras- 
tically change the instability patterns by selecting the 
strong instabilities (large detachments) over the weak in- 
stabilities (small detachments). This is exactly what we 
observe in our numerical studies. 

By applying the general stability analysis of nonlinear 
discrete systemsEZI we obtain three possibilities for mul- 
tiplier collisions. The first is realized when the collision 
takes place on the positive real axis in the complex plane. 
This implies that a multiplier is colliding with its com- 
plex conjugate partner. With (|3^) it follows that 
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Lo = mfi , (33) 

for any integer number m. These are primary reso- 
nances of tlie breather frequency or its higher harmonics 
with any of the frequencies of the EEs. 

The collision of an eigenvalue with its complex conju- 
gated partner can also take place on the negative real 
axis in the complex plane: 

uj = {^^m)VL. (34) 

These are parametric resonances of the breather state 
with EEs . 

The third case is realized when the collision takes place 
away from the real axis. Then a multiplier has to collide 
with a different one, but not with its own complex con- 
jugate partner. It follows that 

wi ± Ci;2 = . (35) 

This is a combination resonance, as the breather fre- 
quency (or its multiple) has to match a sum (difference) 
of the frequencies of two different EE frequencies. 

We stress again that EEs may be localized or delocal- 
ized. While we may proceed with an analytical predic- 
tion of collisions using delocalized Floquet multipliers, we 
have to use numerical calculations to observe collisions 
involving localized ones. Especially the combination res- 
onance may involve either two delocalized, two localized 
or one delocalized and one localized Floquet eigenvalues. 
In addition the above mentioned dissipation-induced se- 
lection of weak and strong instabilities will result in some 
possible collisions being harmless (leaving the breather 
stable) while others will turn out to be important for 
understanding breather instabilities. 

VI. NUMERICAL SIMULATIONS OF 
BREATHER DYNAMICS 

To study the breather dynamics, we performed direct 
numerical simulations of the set of equations (||). All 
simulations were carried out for JJLs with iV = 10 cells. 
We impose open boundary conditions and use the 4th 
order Runge-Kutta method. Time is measured in dimen- 
sionless time units. The initial value of the dc bias was 
7 = 0.8. We choose proper initial conditions that lead 
to the relaxation of the system into a particular breather 
state of left-right symmetry with one resistive vertical 
junction, as in Fig.^ . After a waiting time of 500 time 
units we use the next 500 time units to calculate the time 
averaged characteristics of the state. We then decrease 
the dc bias 7 by a tiny step of A7 = 0.0005 and repeat 
the procedure. We checked that our results do not change 
upon further increase of the waiting time. We varied the 
anisotropy r\ and the inductance of the cell fih while the 
dissipation a = 0.1 was fixed. We will comment on hys- 
teresis effects due to additional increasing of the current 
in the conclusion. 



There are three different ways to monitor the simula- 
tions. The first one is the I — V characteristics, that is 
the dependence of the averaged voltage drop across the 
resistive vertical junction on the dc bias. Furthermore we 
obtain the power Pac (|3^) of ac oscillations of the vertical 
junction at the edge of JJLs. Finally we generate time- 
resolved images (movies) of the full dynamical behavior 
of the ladder so that we can visually check whether the 
system still resides in the initially chosen breather state, 
or switches into another state. Our results are presented 
in Figs. p[]Tl|, where each figure consists of two parts. In 
the left hand parts the I — V characteristics are shown 
(solid lines) together with the approximate results from 
section III (dotted lines). The vertically oriented dashed 
lines indicate the band edges of the linear EWs. In the 
right hand part of the figures we show the dependence of 
Pac (solid lines) on the dc bias together with our approx- 
imate analysis from section IV and Appendix A (dashed 
lines) where appropriate. 

We start with the case of small (3l values. For (3l — 0.2 
and 77 = 1.15 (Figj5|) the breather is easily excited, and 
its frequency is located below io+[q). In the Pac plot we 
observe peaks that are due to the resonance of the second 
and third harmonics of the breather with uj+{q). These 
resonances are primary ones as discussed above. Note 
that their presence is barely seen on the I — V curve. 
The series of observed peaks is related to the finite size 
of the system, and therefore to the resonant interaction 
of the breather with a discrete set of cavity modes as was 
discussed at the end of section II. We tested our interpre- 
tation by increasing the size of the system and observed 
the predicted increase in the number of resonance peaks. 
Close to the lowest possible current (around 7 — 0.55) 
we observe a switching to another breather state, which 
however has the same symmetry and spatial structure. 
Note that shortly after this switching (upon further low- 
ering of the current) we lose the breather and the system 
switches to the superconducting ground state. 

For lower values of the anisotropy r] = 0.35 (Fig||) the 
resonances are again not detectable in the I — V curve. 
However, by monitoring Pac we observe the singularities 
that correspond to the primary resonance 2Q = i-Uj^(q) . 
Moreover, at the dc bias 7 ~ 0.35 we detect a weak third- 
order primary resonance 3Q — uj+{q) in the breather tail. 
The dashed line in the right part of the figure is the pre- 
diction of Pac using our approximate tail analysis. Note 
that our approximate tail analysis is based on the as- 
sumption of a dense spectrum of EWs. Consequently 
the calculated Pac presents an envelope of the numeri- 
cally observed series of discrete peaks. 
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FIG. 5. I — V characteristics and the Pac dependence on 
the dc bias current for a = 0.1, /3l — 0.2, 77 = 1.15. 




FIG. 7. I ~ V characteristics and the Pac dependence on 
the dc bias current for a — 0.1, /3l = 0.5, rj = 1.15. 
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FIG. 6. I — V characteristics and the Pac dependence on 
the dc bias current for a = 0.1, /3l ~ 0.2, 77 = 0.35. 

Next we increased the inductance of the cell to (3l — 
0.5. In Fig.0 we show the results for 77 = 1.15. For 
the initial value of the dc bias 7 = 0.8 the breather 
frequency is already located inside the Lo+{q) band of 
EWs, and this primary resonance is observed in the I~V 
curve. Indeed, the slope of the I — V curve is larger 
than the prediction ( pl|) which does not take into account 
resonant interactions with EEs. With decreasing dc bias, 
the breather frequency is lowered and the above primary 
resonance disappears. However at lower current values 
the next primary resonances 2f2 — LiJ^{q) occur and are 
observable, both in the breather tail and in the I — V 
characteristics. 

The primary resonance structures {fl = a;+(g)) in the 
large current domain are also observed for smaller values 
of the anisotropy parameter rj = 0.5 (Fig.^. In this case, 
they manifesj_.themselves through resonant steps in the 
I — V curve.ta At lower values of the dc bias we again 
observe primary resonances with m = 2. 
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FIG. 8. I — V characteristics and the Pa, 
the dc bias current for a — 0.1, /3l = 0.5, rj - 
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For ~ 1-0 and rj = 0.5 (FigJ|) the breather fre- 
quency is located above uj+iq) for large current values. 
Upon decreasing the dc bias we observe a peculiar switch 
to a different breather state with the same spatial struc- 
ture but a lower frequency located inside u}+{q). The 
most interesting feature here is that shortly before the 
switching the breather frequency is clearly larger and 
outside of the uj+{q) region. Upon further lowering of 
the dc bias we observe primary resonances = Lu^{q), 
and corresponding resonant steps m I ~ V curve. 
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FIG. 9. I — V characteristics and the Pac dependence on 
the dc bias current for q = 0.1, /3l = 1.0, r\ = 0.5. 



Let us increase the inductance of the cell /3i even fur- 
ther. For = 3 and rj — 0.5 (Fig.|o|) we again find 
that breather frequencies are located above uj+ for large 
current values. Similar to the previous case we observe 
a switching when the breather frequency is clearly out- 
side (above) the branch uj^(q). This brings the system 
into another breather state with the same spatial struc- 
ture, but with a frequency again located above u;+. This 
highly nonlinear state is then lost by switching to the su- 
perconducting ground state after further decrease of the 
dc bias. 
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FIG. 10. I ~ V characteristics and the Pac dependence on 
the dc bias current for a — 0.1, I3l ~ 3.0, rj = 0.5. 

Now we come to an interesting observation. Lowering 
the anisotropy rj = 0.35 (Fig.[ll|) we again observe the 
switching at a breather frequency being located above 
However the switching increases the voltage drop. 
The new state is of different internal structure. We re- 
mind the reader that all previous numerical results have 



been obtained for a breather with a structure as in Fig.gb. 
Here we find that after the switch the new breather state 
is characterized by three vertical junctions being in the 
resistive state. At the same time the symmetry is broken. 
In fact the new state exactly corresponds to the example 
given in Fig.||d. Note that similar switchings (which lead 
to an increase of the number of resistive juuctipns) have 
been reported in early experimental studies-cHlj The left- 
right symmetry breaking leads to new interesting features 
in the breather tails. The no-symmetry breather has two 
times lower frequency than the voltage drop across the 
vertical junction. Thus the new breather frequency is in- 
side the upper band (jj+{q), and a primary resonance is 
clearly observed in the Pac dependence on 7. 

In order to test the influence of small fluctuations on 
the I —V curves we repeated the simulations in the pres- 
ence of small noise with amplitude w 10^^. All obtained 
results are stable except the switching outcome in Fig.pT| 
. While this switching occurs at the same current value, 
the new breather state is changed. In particular we ob- 
served the left-right symmetry breather with three junc- 
tions being in the resistive state. Thus we find extreme 
sensitiveness of the outcoming breather structure (includ- 
ing its symmetries) to small fiuctuations. This implies 
that the boundaries of the volumes of attraction of differ- 
ent (breather) attractors are entangled in a very peculiar 
way. We can reliably predict the switching position, but 
not the outcome of the switching. 

We also numerically simulated the breather dynamics 
in the JJL with an extremely large inductance of the cell 
= 500. We did not find any indication of resonances 
and instabilities. We argue that the reason for that is the 
weak dispersion of the linear EWs for such large values 
of (3l- This implies that interactions along the ladder are 
weak. The breather is continued to small current values 
until it switches to the superconducting state at the dc 
dc bias 7 = 0.22 for rj — 0.35. This particular value can 
be obtained by making use of the simple dc analysis ( |T^ ) 
and the standard theory of thfi-jretxapping current in a 
single small Josephson junctiont2Hl3 



7r = (1 + 2r,) 



4a 



(36) 



This equation yields a value of 0.22 for the considered 
case, in good agreement with the numerical observation. 
Note that within this theory retrapping occurs purely 
due to energy considerations, not due to resonances (or 
instabilities). 

It is very important to notice that for cases with small 
or intermediate values of inductance of the cell (Fig.||- 
[Tl] ) the observed currents at which we lose the breather 
state and switch to the superconducting one exceed the 
expected retrapping values (|3^). We will explain this 
disagreement in the next section. 

Motivated by the above findings we investigated the 
loss of the HWS upon lowering the current. We recall 
that in this state all vertical junctions are resistive and 
all horizontal ones are superconducting. Usually it is 
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assumed that the HWS loss is again due to a standard 
retrapping mechanism. It is important that any numeri- 
cal simulation of such a process is done with the addition 
of some weak noise, because the processor will otherwise 
perform a perfect simulation of a single junction repeated 
N + 1 times. We chose = 3 and 77 = 0.35. The I-V 
characteristic is shown by a thick dashed line in Fig.|ll| . 
The expected retrapping current Aa/n = 0.127 is clearly 
not reached. Instead we observe the loss of the HWS at 
7 = 0.273. At the same time it follows from equations (H) 
that the HWS exists as a solution down to the retrapping 
current of a single junction, i.e. down to 7 = 0.127 ! So 
in this case we conclude that the numerically observed 
loss of the HWS at 7 = 0.273 is due to an instability. 
The HWS continues to exist as a solution down to the 
standard retrapping current, but it is an unstable state. 
This novel result is very important, since very often the 
current value of HWS loss in the absence of a magnetic 
field is used to estimate different parameters of the sys- 
tem assuming that the HWS is behaving similar to a 
single junction. Our results show that this is definitely 
not the case. 
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FIG. 11. I — V characteristics and the Pac dependence on 
the dc bias current for a — 0.1, /3l — 3.0, rj = 0.35. 



VII. EVALUATION OF RESONANCES AND EWS 
SPECTROSCOPY 

This section is devoted to a quantitative explanation 
of the observed resonances and switchings. 



A. Primary resonances 

Primary resonances are characterized by mfl ~ lu 
where lu is some EE frequency. We detected various pri- 
mary resonances with extended EWs. The case m — I, 
which corresponds to the breather frequency being lo- 
cated inside the uj^(q) band, shows up with resonant 



steps in the I — V curves (see. Figs. |^^. The finite 
number of observed resonant steps is due to the discrete 
spectrum of the excited cavity modes. In addition we 
observe strong variations of the breather tail amplitudes. 

Higher order primary resonances (m — 2, 3) are much 
less pronounced in the I—V characteristics. They mainly 
lead to a weaker localization of the breather tail and can 
be clearly detected in the form of sharp peaks in the 
Pacil) dependence. Since the breather in our case has 
left-right symmetry, the only linear cavity modes which 
can be exited are symmetric (see section II) . These modes 
are characterized by even values of k (p^). We start our 
evaluation of these resonances with the case shown in 
Fig.^. We determine the 7 value of each observed peak 
in Pacil) and thus obtain the corresponding breather 
frequency ^2(7). We then compare its multiples with 
the discrete spectrum of linear mode frequencies of the 
Lu^ branch. The numbers are listed in Table I. We find 
that all observed resonances are due to symmetric linear 
modes (even fc = 2, 4, 6, 8, 10) as expected. 





k 




20, 


30 


1 


5.475 






2 


5.581 


5.594 


5.590 


3 


5.745 






4 


5.949 


5.947 


5.945 


5 


6.174 






6 


6.400 


6.387 


6.384 


7 


6.609 






8 


6.788 


6.769 


6.772 


9 


6.923 






10 


7.008 


6.960 


6.958 





TABLE I. Comparison of theoretical predictions of pri- 
mary resonances with the numerical results from Fig^. First 
column: k. Second column: the spectrum of EWs u!+{qk) (cf. 
(p^. Third column: 2f2 obtained from the peak positions 
of Pac (7) in Fig.^(a). Fourth column: 3f2 obtained from the 
peak positions of Pac(7) in Fig.^(b). 
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The same method of analysis allows us to conclude 
that the three peaks in Fig.|| (in decreasing order of 
dc bias) are due to the following resonances: (m, k) = 
(2, 6); (2, 4); (3, 8). Similarly the shoulder and the peak 
in Fig.||are due to resonances with {m,k) — (2, 10); (2, 8). 
Finally the resonances in Fig.^ correspond to the values 
(m, k) = (2, 8); (2, 6); (2, 4). Note that in all of these cases 
the deviations between the theoretical and observed num- 
bers are much less than the frequency difference between 
adjacent cavity modes. 

B. Parametric resonances 

So far we did not comment on the nature of the switch- 
ing from a breather state to the superconducting state 
for small and intermediate values of P^. The Floquet 
analysis results show that all these switchings are due 
to an instability of the breather. In terms of Floquet 
multipliers, all of these instabilities are due to a collision 
of two localized Floquet multipliers on the negative real 
axis. The breather state continues to exist as a solution 
to the dynamical equations for lower current values, but 
it is unstable. Note that the so-called retrapping mech- 
anism instead (as for a single junction) uses the critical 
current value as a criterion for retrapping. This argu- 
ment is based purely on energy considerations and does 
not take into account any resonance mechanism. This is 
not surprising, as a single junction has no other degrees 
of freedom it may resonate with. Below the retrapping 
current, the resistive state disappears in this case. So 
we may state that the switching from a breather state 
to the superconducting one as observed in our simula- 
tions is usually driven by resonances with localized EEs 
(frequency matching) and is not due to energy effects 
(current value matching). 

C. Combination resonances 

Let us discuss the nature of the switchings of the 
breather for intermediate Pl values when the breather 
frequency is located above the branch Lu^{q). These 
switchings are again due to an instability. It is char- 
acterized by Floquet multipliers colliding away from the 
real axis. As discussed in Section V this corresponds to 
a combination resonance. The numerical Floquet analy- 
sis shows that one of the two participating multipliers is 
a localized one (which bifurcates from the lower branch 
W- (<;))., while the second one belongs to the delocalized 
spectrum of uj+ (q) . The Floquet multiplier which finally 
leaves the unit circle is a localized one. So again the 
instability of the breather is driven by a localized pertur- 
bation. 




FIG. 12. Arguments and absolute values of Floquet multi- 
pliers versus 7 for the breather state in Fig.pl] 

In Fig.^ the dependence of the arguments and abso- 
lute values of all relevant Floquet multipliers is shown for 
the breather of Figjl^. For convenience we do not plot 
the complex conjugate multipliers and restrict the argu- 
ments to < arg(fi) < n. The narrow band uj-{q) and 
broad band uj+{q) are nicely observed. The degenerate 
band ujq is located slightly above uj^ (q) . This band does 
not interact with other multipliers when crossing them, 
as expected from our analytical considerations. The two 
separated arguments which are located below the i-u^{q) 
band have localized eigenvectors .Ell 

In the plot of the absolute values we observe the pre- 
dicted values n = 1 and e"""^*". The multipliers which 
correspond to lines between these two states generally 
reside on the circle with radius (pi]). Many of them de- 
part from this circle due to collisions. At current val- 
ues of 0.55 < 7 < 0.7 we observe parametric resonances 
2il = uj^{q), which belong to the set of weak resonances 
and do not evolve into a global instability. However it 
is possible that a slight variation of control parameters 
(e.g. decreasing the damping a) might change these res- 
onances into strong ones. Then we would expect sudden 



11 



instabilities of the breather state at these large current 
values. In our case the global instability is realized when 
one of the localized multipliers collides with the oj+{q) 
band around 7 — 0.48. Subsequent lowering of the cur- 
rent leads to a fast escape of this multiplier from the unit 
circle and to the observed switching. 

The importance of localized EEs for the destabiliza- 
tion of a breather is simply due to the localized nature of 
the latter. It is hard (if not impossible) for a breather to 
generate a parametric instability through extended EWs 
alone, as these excitations are damped out far from the 
breather center. In contrast localized EEs do not travel 
away from the breather center. These modes can be ef- 
fectively excited by the breather, leading to an instability 
of the latter. 

To understand the nature of the observed instability 
of the HWS we show a similar Floquet multiplier plot in 
Fig-H. 
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FIG. 13. Arguments and absolute values of Floquet multi- 
pliers versus 7 for the HWS in Fig.pj|. 

As in the previous plot, we observe weak paramet- 



ric resonances of the upper EW band at current values 
7 ^ 0.38 which do not evolve into a global instability. 
Again the observed instability is driven by a combination 
resonance at 7 = 0.175. Since the HWS is an extended 
state, all Floquet multipliers are also extended. The com- 
bination resonance is due to the collision of two Floquet 
multipliers belonging to the two EW branches oj± {q) with 
q = IT (we remind the reader that the EW spectrum of 
the HWS is different from that of the superconducting 
state and can be obtained by putting 7 = 1 in (p^). 
Note that indeed for the present case, the frequency of 
the HWS at the instability equals 2.73, while the value of 
the combination i-U-(tt) + w+(7r) = 2.69. For current val- 
ues 7 ^ 0.18 the HWS is becoming stable again. However 
around 7 — 0.17 another even stronger instability due to 
parametric resonance sets in, which brings the HWS to 
the next instability well above the expected retrapping 
current. 



VIII. CONCLUSION 

We have presented analytical and numerical studies of 
breather properties in Josephson junction ladders. Our 
results confirm and substantially extend early sugges- 
tions that breathers may resonate in different ways with 
localized and extended electromagnetic modes. The nu- 
merical studies have been done in a parallel manner to 
the way experiments are conducted. The variation of the 
control parameter 7 allows one to continuously change 
the breather frequency, whereas the linear mode spec- 
trum is not significantly changed in the domain of inter- 
est. We observed primary resonances mfl = uj+{qk) with 
extended EWs, parametric resonances Vl — 2uj with local- 
ized EEs, and combination resonances mD, = lo -\- LOj^{qi) 
with a localized EE and a delocalized EW participating. 
We also observed a combination resonance which leads to 
a switching from a small breather (one resistive vertical 
junction) to a larger one (three resistive vertical junc- 
tions) together with a possible symmetry lowering of the 
breather. 

The primary resonances with extended EWs lead to 
singularities in the breather tails. This allows one to de- 
velop a spectroscopy of EWs by monitoring Pac versus 7. 
Such a spectroscopy can be experimentally realized, e.g. 
with the help r©f|-a, well known Josephson junction detec- 
tor technique .I12IE3'E2I It could be important for obtaining 
a coherent source of high frequency radiation, since in 
such a resonance the whole breather tail starts to co- 
herently oscillate with large amplitudes. The resistive 
breather center serves as a region of energy input via a 
dc bias. 

Our studies show that the main control parameter (in 
addition to the dc bias) is the self inductance (3l- For 
small Pl values the breather frequency is located between 
the two branches of EWs, LO±{q). One may perform spec- 
troscopy of EWs of the upper branch, or observe para- 
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metric instability of a breather due to localized modes. 
Moreover as increases, resonant steps in the I — V 
characteristics can be observed. For intermediate (3l > I 
values, the breather frequency is located above uj^{q). 
In this case we observe combination resonances due to 
localized and delocalized EEs, which may result in a un- 
usual sensitivity of the switching outcome on small fluc- 
tuations. Large Pl values stabilize the breather states, 
make resonances impossible and lead to a standard re- 
trapping mechanism for breather switching to the super- 
conducting ground state. This is likely the situation for 
the reported experimental data in [Rcf. ^ . However lower 
(3l values allow for the appearance of the above listed 
resonances (and perhaps even other still unobserved res- 
onances). We believe that our findings will help to make 
the proper parameter choice when designing new ladders 
for experiments. 

Note that throughout our studies we always decrease 
the dc bias 7. Let us consider a breather state Nr.l 
which becomes unstable upon lowering the dc bias at a 
certain value 71 . Let this be a case where the system will 
switch to another stable breather state Nr. 2. This new 
breather state is in fact keeping its stability not only upon 
further lowering of the dc bias, but also upon a reversing 
(increasing) of the dc bias. Thus we find that there ex- 
ist dc bias windows in which both the starting breather 
state Nr.l and the new breather state Nr. 2 are stable. 
Even though their dc spatial structures (cf. Fig.||) may 
be identical, the average voltage drops (and frequencies) 
are in general different. Further increasing of the dc bias 
while staying on the new breather state Nr. 2 will lead 
to an instability and switching at 72 > 71. In case the 
switching brings the system back to the breather state 
Nr.l we are faced with the well-known two-state hystere- 
sis phenomenon in Josephson junction systems. However 
we also observed cases when the switching due to an in- 
stability of breather state Nr. 2 (upon increasing the dc 
bias) brings the system to yet another state, which dif- 
fers from breather state Nr.l, e.g. simply to the HWS. 
In such a case the hysteretic behavior is of more complex 
nature. To keep the discussion of our studies as clear as 
possible we did not present data for increasing current. 

Another observation is that the expected vakies of the 
retrapping current based on a pure dc analysisll3 are too 
low to match the observed values at which the breather 
switches to the superconducting ground state. Only for 
very large Pl values do we observe agreement. For all 
other cases the breather switches to the superconducting 
state via an instability driven by parametric resonance 
with localized EEs. Moreover the HWS also undergoes an 
instability which is due to combination resonances with 
extended EWs. 

In this work we always started with a breather con- 
figuration as in Fig.||D, at large current values. It may 
be expected that the results for other starting configura- 
tions will show up with similar properties, and perhaps 
with additional types of resonances as well. This may be 
due to the fact that the structure of the phase space is 



very complex, being separated in many different regions 
of attraction of different attractors. It is this complex- 
ity which makes the understanding of breather properties 
both a fascinating and complicated undertaking. 
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APPENDIX A: POWER OF AC LIBRATIONS AT 
THE EDGE OF A JJL 

Here we derive the time-average power of ac librations 
at the edge of a JJL. This characteristics is proportional 
to the average kinetic energy (0^/2). In order to obtain 
an expression for the kinetic energy, we have to determine 
the dynamics of the junction at the JJL edge. For this 
we write the system of equations (p5|) in a matrix form: 



Av = . 

Here, v is an unknown vector 



Ah 



and A is a 2 X 2 matrix 



i(A) 



aii(A) ai2(A) 
021 (A) 022 (A) 



(Al) 



(A2) 



(A3) 



Pi 



■ cosh A, ai2(A) 







Ml 



B. A nonzero solu- 



where an (A) 

a2i(A) = -■gi^(l-e^), and a22(A) 

tion exists if the determinant of A vanishes: det A(Xo) = 
0. The parameter Aq was determined in section IV ([2^). 
The components of the vector v satisfy the condition: 



022 
a2i 



(A4) 



To determine the components V and H separately we 
have to impose an additional condition at the breather 
center. This condition is not known exactly due to the 
complex dynamics in the resistive breather center. Nev- 
ertheless the ac librations in the breather tails are weakly 
depending on it. Here we use the simplest normalization 
condition 



|AJ2 + |AJ2 = 1. 



(A5) 



Substituting (A4) into (A5) we obtain 
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^21 



V|a2iP + |a22P 

022021 



V|a2i|* + |a22P|a2iP 



(A6) 



Due to the up-down symmetry in the breather tail, the 
dynamics of Josephson phases at the edge of the JJL can 
be written in the form 



(a22a2i)e^^^°^" cos(a(Ao)n + nt) 
v/|a2i|* + |a22P|a2iP 
3(a22a2i)e**(^'')" sin(3(Ao)n + m) 
^ v/|«2i|4 + |a22PR^ 

. |a2ipe'^(^'')"cos(Q(Ao)n + m) 

<^n = — 



(A7) 



V|a2ip + |a22p 

We finally obtain the expression for the average kinetic 
energy (for n < 0) 



1 



^ |a22a2i| e 



2„Si(Ao)2n 



4(|a2i|4 + |a22P|a2iP) 



(A8) 



APPENDIX B: LINEAR STABILITY OF THE 
BREATHER IN JJL. 

The stability of periodic motion is analyzed with the 
help of the Floquet theoryO'cil Linearizing the system 
around a time periodic breather solution, we obtain 

el + ael + Al{t)el = ^{Ael + Ve^i - Ve^i) 

PL 

■£l + ae^^ + A':Si)e^^ = -^(Ve'^; + e,'j - e^) (Bl) 

VPL 

ll + all + Ai{t)ll = ^(Ve^; + - et) , 

where An it) are time-periodic coefficients determined by 
the given breather state. 
The substitution 



(B2) 



allows one to eliminate the dissipation 



1 



< + = — (A< + V«^i - Vk^i) 

PL 



(B4) 

where Bn{t) = —\ct^ + An{t). Introducing new variables 



we find the system of equations 



(B5) 



PL PL^A? 

z^: -I- = -^(^n - 5:;) - ^v^;; (b6) 

?7PL Plv^ 



/3lV^ " 

These equations describe a Hamiltonian system, namely 



dpn 

dn 

dZn 



(B7) 



where z„ = , Pn = and the 

Hamiltonian 7i(0„,p„,t) is 

n 

n 



(B8) 



2/3l?7 



Since the particular Hamiltonian can be represented in a 
general quadratic form, the symplectic product of two dif- 
ferent trajectories ^^^(i), z„(t)} and {Pni^^) , z^{t)} does 
not change in timetll 



^ = Y.iPnit)Zn{t)-Pn{t)z;{t)] . 
n 

Rewriting our set of equations in the form 

0Z„. = -T^Opn + „^ 0Z„ 



(B9) 



-Sp„ 



dp^ dpndZn 

-OPn + ^T^OZ„ 



ZnOPn OZn 

and using the notation 

^ \-E 
where E is the identity matrix, we obtain 



^Pn \ _ j-w^n ^ 



5z 



6Zn 



(BIO) 



(Bll) 



(B12) 



14 



where V^Ti is the Hessian of H. 

Let us c onsid er the following map by integrating the 
equations ( B12| ) over one period of the initial solution 



Spn{Tb) 



uin) 



SpniO) 



(B13) 



Since the form J is symplectic (B9) , U{Ti,) is symplectic 
too. As a result, we find that the eigenvalues of U{Ti,) 
have to fulfill the condition that if i/ is an eigenvalue 
then V* and 4r are also eigenvalues. Note that for a 
marginally stable periodic motion of a Hamiltonian sys- 
tem the Floquet eigenvalues v are located on the unit 
circle. Switching to an unstable state is realized by colli- 
sions of eigenvalues on the unit circle and departing from 
it. 

All of the obtained relations for eigenvalues (and eigen- 
vectors) can be rewritten for the original 5n{t) variables 
(call the corresponding Floquet eigenvalues /i). First of 
all, we expect that most of the eigenvalues will be located 
on a circle of radius (^Tj) . Further it follows that if ji is an 



eigenvalue then 



and 



are eigenvalues too. 



So, all scenarios of collisions of the eigenvalues are simi- 
lar to the Hamiltonian case. For a periodic motion to be 
stable we need all eigenvalues /i to have absolute values 
less or equal to one, i.e. the complex numbers should 
reside inside the unit circle in the complex plane. Since 
the radius R{a) < 1 it needs further finite variation of the 
control parameters (after a collision) to enforce an eigen- 
value to traverse the finite distance to the unit circle and 
to exit it. Thus not every collision will lead to an instabil- 
ity. However, the decrease of the dissipation parameter 
a may tune the system closer to the Hamiltonian case, 
making its states more sensitive to any occurring colli- 
sions. Note that due to the periodicity of the breather 
state there is always one eigenvalue ^ — 1, whose eigen- 
vector is tangent to the breather orbit. Consequently 
we always find another eigenvalue /x = e"""^' which is 
located on the positive real axis inside the unit circle. 
During all numerical computations of Floquet eigenval- 
ues and eigenvectors the above properties were tested and 
complete agreement was found. 
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